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Abstract 

High energy factorization predictions for F2 are derived using BFKL 
descriptions of the proton structure function F2 at HERA. The model 
parameters are fixed by a fit of F2 at small x. Two different approaches 
of the non perturbative proton input are shown to correspond to the 
factorization at the gluon or quark level, respectively. The predic- 
tions for F2 are in agreement with the data within the present error 
bars. However, the photon wave- function formulation (factorization 
at quark level) predicts significantly higher F2 than both gluon fac- 
torization and a next-leading order DGLAP model. 



1 Motivation 



High energy factorization |]T], is a QCD factorization scheme suited for high- 
energy hard processes - and in particular for deep-inelastic e'^—p scattering 
in the small x-regime (x ~ ^) {W^ is the energy squared in the virtual 
photon-proton center of mass frame). This scheme takes into account the 
resummation of the (^a^log^^ terms in the QCD perturbative expansion 
of the structure functions. It amounts to proving the factorization of the 
perturbative amplitude in terms of a gluon-gluon BFKL kernel P| convoluted 
with a first order (virtual)photon-(virtual)gluon cross section. A graphical 
illustration of this property is given in figure 1. The main difference between 
this scheme and the renormalization group factorization which is valid at 
finite x is that the former involves a bidimensional integration in transverse 
momentum k±, whereas the latter is a convolution in energy. Interestingly 
enough, high energy factorization can be applied to the resummation of the 
(^ag log terms in heavy quark pair photoproduction, and the combined 

resummation for heavy quark leptoproduction involving lag log -^1 and 

(as log terms can also be studied [||]. Our aim is to test the high energy 
factorization predictions for the massive charm quark pair contribution to F2 
at HERA. 

In references ||, a formulation and a phenomenological analysis of the 
proton structure function F2 satisfying high energy factorization in the frame- 
work of Mueller's color dipole picture [0 of BFKL dynamics were proposed. 
F2 is the triple transverse momentum convolution of a coefficient function, 
the BFKL kernel and a non perturbative term. The non perturbative inputs, 
which can be interpreted as the density and size of primary dipoles in the 
parton sea of the proton 0, lead to a satisfactory description of F2 at HERA 
range, with 3 parameters. In this analysis, quarks were assumed massless, 
and the convolution integral was approximated by a steepest-descent method 
for sake of simplicity. However, since the coefficient functions were consid- 
ered constant and the quarks massless, the procedure was not specifically a 
check of high energy factorization itself. 

In the following, we shall discuss the high energy factorization predictions 
by taking full account of the convolution integral and of the massive quark 
component in the phenomenological analysis. 

On the experimental point of view, the HERA experiments have recently 
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published data for the contribution of charmed meson production to the 
structure function F2. Their analysis was based on and D* meson tagging. 
This allows to single out the charm contribution F2 to the total structure 
function and thus to investigate the quark mass dependence of the struc- 
ture functions. We use this nice possibility to check whether the high energy 
factorization scheme gives a correct prediction for the mass-dependent contri- 
bution both in X and Q^. Our analysis can be considered as the analog within 
the BFKL dynamics of the one which can be performed using ordinary renor- 
malization group evolution |P, |10|. Thus, we are also aiming to see whether 



there are agreements or differences in charm leptoproduction between the 
different schemes available for analyzing deep-inelastic scattering. 

However, the high energy factorization scheme adapted to the proton 
structure functions is not uniquely defined in terms of the separation be- 
tween perturbative and non-perturbative contributions. Indeed, an a priori 
different factorization exists in which the virtual photon is described in terms 
of a quark- ant iquark wave function configuration which then interacts with 



the proton (see for instance |TT] , |T2[). One of our goals is to investigate the 
possible differences between the two schemes. Recently also, some theoreti- 
cal doubt has been cast on the validity of the operator product expansion in 
the low-x range due to the fc^-diffusion property [jl3l inherent to the BFKL 
dynamics. In the context of our structure function studies, we thus also 
want to discuss phenomenologically the kinematical range in x and where 
high energy factorization is valid and not spoiled by /c^-diffusion in the low 
momentum region. 

The plan of our study is the following. Section 2 recalls the formulation 
of high energy factorization for the pair production of quarks of mass M in 
terms of the coefficient functions and unintegrated gluon distribution in the 
proton. We give the M-dependent expression of the coefficient functions. In 
section 3, we derive the expression of the proton structure function including 
the charm component in the framework of high energy factorization. We de- 
rive a new constraint on the unintegrated gluon distribution function due to 
the renormalization group property at high k^. In section 4, we consider the 
alternative model of factorization based on the wave function formulation. 
We show the connection between the two models. We prove the equivalence 
at perturbative level between both schemes for the hard vertex while they 
are based on different non-perturbative inputs. In section 5, both models are 
applied to a phenomenological fit of F2 and parameter free predictions for 
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F2, Fl, and F^. The final section 6 contains our conclusions on the validity 
of the high energy factorization predictions at HERA and an outlook on fu- 
ture studies. Taking into account both present experimental and theoretical 
uncertainties, high energy factorization is proved to be in agreement with the 
present experimental data. However, more precise data could discriminate 
between the different schemes, since the photon wave function formulation 
predicts higher than both the gluon factorization and a next-leading order 
DGLAP model. 



2 High energy factorization 

Let us compute the pair production of quarks of arbitrary mass M at the vir- 
tual photon vertex of deep-inelastic scattering (see figure 1). Using the high 
energy factorization scheme the inclusive transverse (resp. longitudinal) 
structure functions Ft (Fl) can be expressed as follows: 

ir,.T.L(y-V, 11/M, kJM) T(y, fci) , (1) 

where = —q^ is the virtuality of the photon. Y represents the rapidity 
range available for the reaction. J-'{y, k±) is the unintegrated gluon distribu- 
tion m, which describes the probability of finding a gluon with a longitudinal 
momentum fraction z = and two-dimensional transverse momentum k± 
in the target. q± is the photon transverse momentum in the photon-proton 
centre of mass frame. d'^*g^T,L is the hard cross section for (virtual)photon- 
(virtual)gluon fusion computed at order agaem- 

In order to express more easily high energy factorization, it is convenient 
to work with triple Mellin-transforms with respect to the rapidity Y and 
transverse momenta k± and q±. The result is: 



Ft,l{Y,Q\M^ 




d'y [AM 



2 \ 7 



2i7r7 V g2 



/ / 

2^ 



iv(7',7)-^iv(7;Q?), (2) 



where stands for the factorization scale and the inverse Mellin-transforms 
in 7, 7' and are given by complex integrals along the axes | ± ioo, \ ± ioo 
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and No ± ioo {Nq larger than the real part of the rightmost singularity in 
the iV-plane) respectively. By definition, hT,L,N{'j' ,j) is the triple inverse 
Mellin-transform of the hard 7*-gluon cross-section, namely: 



^n{1] Q}) is obtained from the unintegrated gluon distribution J-'{y, k±). It 
reads: 

^„(0.Qj).f #e-.-/^(|)'Vfe.*J. (4) 

Note that the coefficient functions /iT,L,Af are known to have a weak depen- 
dence on III, hence we will consider in the following only their values at 
= denoted /it,l(7; M^). 

The final expression for the high energy factorized structure function is 
FT..(r,Q,M) = y — (^^j /^T..(7;M) (5) 

where 

Hy,r.Q)) = J^^'^''^^(r,Q}) (6) 



and the coefficient functions related to Fx and respectively are [|T4 
as 4-^ r(l+7)r3(l-7) 1 



67r(l + 27)(l- |7) r(2-27) 1 + 



4Af2 



(37-1) + -^(7 -2) + 1 + 



2(l+7)(2-7) + -^(7+7-67'— ^(7-2))j x 



1 3 

2' 2' 1 + 
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and 



hLir,M' 



2 ^ c^s 4-^ r(i+7)r^(i-7) 1 

67r(l + 27)(l-|7) r(2-27) 1 + ^ 



2(1-7) + 3-^)+ 1 + -^ X 



U(l-7)-^(3^ + 4(l-oO) X 



1 3 

2' 2' 1 + 



2^1 1 1-7,;^;;^;^^ 1 ) , (8) 



where 2-^1(0, h\ c; () is the hypergeometric function |Tf 

It is easy to check that in both hmits (hght quarks [|T^) and 

Q"^ ^ (heavy flavour photoproduction [|l|), the well-known expressions for 
the coefficient functions are recovered. 

Inserting the expressions (7,8) for /iT,L(7; M^) and a model for J-'{Y, 7; Q'j) 
in formula (^, and summing over all active flavours weighted by their electric 
charges squared, we obtain the explicit expression for the proton structure 
functions. 



3 The proton structure functions 

Let us now introduce the QCD dipole model for the proton structure func- 
tions. This model satisfies high energy factorization and gives a phenomeno- 
logical description of the unintegrated structure function J^n{.1] Q}), see (4), 
in the framework of BFKL dynamics. We parametrize the factor 

•^7v(7;Q/) ^ ^7v(7; Q}) . . 

7 N-^xil)' 

Notice the factor I/7 which corresponds to the k± integration in Mellin 
transform and thus, u!]\f{'y;Q'j) appears as the residue of the BFKL pole in 
the integrated gluon distribution in the target. J^n is assumed to contain the 
well-known BFKL sing ularity [| at iV = ^^(7) with 

X(7) = 2^{1) - M/(7) - ^(1^) . (10) 
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The function c<j7v(7;Q/) will explicitly depend on the nature of the target. 
It may also contain other singularities than the BFKL pole, which will be 
discussed later on. For a target of small size (parton, massive onium), it can 
be deduced from perturbative calculations, whereas for extended targets like 
a proton, it is supphed by a model. 

Following the suggestion of ref. 0, one assumes for the proton structure 
functions the following scaling form: 

^7v.o(7;Q?)=^(7)(^) ' , (11) 

where uj{'y) and Qq are the non-perturbative inputs of the model. Indeed, 
the scaling assumption (|TT|) is the simplest one allowing to obtain formulae 
independent of the arbitrary factorization scale Qf (see equation (||)). Note 
that Qo can be considered to be independent of the quark mass M, which 
is quite a reasonable assumption for a non perturbative proton parameter. 
Inserting (P) and the expressions for the coefficient functions (7,8) in equation 
(^, and taking into account the scaling (|TT|), the overall formula reads: 

FtA^, Ql) = ^ el FtAYq, Q\ Ml- Ql) 



Q 



^(7) 



where 



2t7r \QIJ 2tn N - ^xil) '(12) 



^2 = log--^^ (13) 



is the maximal available rapidity range for the produced gluons in associa- 
tion with quarks of mass Mq. The summation J2q takes into account the 
contributions of all active flavours with charge cq, mass Mq. 

The unintegrated gluon distribution in the proton J-'n{i', Q}) is thus 
model dependent in particular through the input function uj{'y). However, let 
us show that it obeys a theoretical constraint when 7 goes to infinity fl^ . In- 



deed, as we shall demonstrate, this limit corresponds to a situation in which 
the intermediate gluon emitted from the dipole (see fig.l) has a large trans- 
verse momentum k±. Hence its evolution from Qq up to fc^. is governed by 
the gluon-gluon DGLAP evolution equation (in its small-a; approximation). 
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Let us then consider formula (4) and single out the integration region 
kj_ ^ . In this hard region where fc^-ordering of intermediate gluons is 
valid, the unintegrated gluon distribution is simply related to Q{Y,k'j_), the 
gluon distribution function in the proton at the scale k^, by: 



Hy, k] 



1 d 



TT dki 



■G{Y,k] 



(14) 



Introducing the Mellin transform in rapidity, the gluon-gluon DGLAP evo- 
lution equation at small x (small N) reads 

as Nc 

dYe''''^g{Y,kl)^g^{kl)^[^^^ g^{Ql). (15) 

Inserting formula (|15D into equation (4) (for k^ ^ Ql), one gets a contribu- 
tion 



f 



\7\~*oo 



1 - 



ttN 



On{qI) , 



(16) 



where 7 = is the well-known DGLAP singularity when N goes to 0. 

As well-known, when 7 is small xil) — I/7 and the BFKL pole (^ can 
be approximated by 



7 



N 



ttN 



N 



(17) 



and thus coincides with the DGLAP pole (|T6|). However, in the large I7I 
region, |x(7)| — 2 log I7I and thus the DGLAP pole at 7 = is separated 
from the BFKL singularity and dominates over it at small N. 

Integrating over this dominant singularity in the inverse Mellin trans- 
form (P), we observe that the unintegrated structure function satisfies the 
constraint 



'Qiy j^i^,Y;Q}) 
'SI 7 



|7|-*oo 



-e ^-y y csN^ {(Jq) 



7r7 



asNc 
7r|7P 



(18) 
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assuming a regular input QN^oiQo)- other hand, the large I7I be- 

haviour of the coefficient functions is given by 

hrAr, M-) ^ ,F,{l-j, 1/2; 3/2; 1) x 



X < terms of the form 7" 



J 

~ e"'^''''' X {Power-like terms} , (19) 

and thus gives an exponential cutoff at 7 ^ ^. 
The following remarks are in order: 

• The obtained behaviour (|r^) for J^n{1'-, Q/) at I7I large is in agreement 

OLsNc 



with the scaling assumption ([Tl]) at the pole 7 



Most importantly, the large I7I behaviour of the coefficient functions 
\hT,L{l I — e"'^!'''! dominates over the unintegrated structure func- 
tion T{Y, 7; Qj) which is only power-like in 7, showing that the main 
integration region is for finite I7I < ^. 

The large I7I behaviour of 7;(5j) is actually not dominated by 
the BFKL singularity but by the DGLAP singularity which differ in 
this domain. However, we shall neglect this modification occurring in 
a domain where the integrand is cutoff by the coefficient functions. A 
more detailed analysis of high energy factorization could need taking 
care of this modification. 



4 High energy factorization at the quark level 

In the previous section, we have found that the main integration region is for 
finite I7I. In fact, the structure of the integrand singularity at 7 = appears 
to be essential both for the theoretical analysis and the phenomenological 
application. We note that in formulae (7) and (8) for the coefficient functions, 
2^1(1--)', |; |; 1/(1+4MVQ^)) oc I/7 when 7 goes to 0. We thus find that the 
Mellin-transform /ir,L,Af(7; ^^)/7 of <5" in formula (1) has in general a double 
pole at 7 = (but for e.g. when = 0). It is easy to realize that this 
double pole corresponds to an extra {logk^/Q"^) in the small fc^-behaviour 
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of the hard photon-gluon cross section due to the quark propagator (see 
fig.l). This behaviour is thus characteristic of the high energy factorization 
formahsm at the perturbative level. 

However, it is well-known that there is an ambiguity in the separation 



between perturbative and non perturbative contributions in the small-k± do- 
main. The relevance of the perturbative double pole depends on the physical 
picture of the non perturbative input. For definiteness, we will consider two 
classes of models relying on different hypotheses on the behaviour of the 
residue function uj{'j) near 7 = 0. In the first type of models, with 

07(7) ~ (constant) , (20) 

7-^0 

one keeps the full perturbative information on the 7*-gluon vertex. This 
corresponds to the factorization at the gluon level (see fig.l). In the second 
class of models, the perturbative singularity at the hard vertex is smoothened 
by the proton scaling function (|TT]), namely 

u;(7) ~ (constant) x 7 . (21) 

7-^0 

The resulting single pole at the hard vertex may be interpreted as a direct 
pointlike coupling of the virtual photon to the quark. This may be interpreted 
as a factorization at the quark level (see fig.l). 

In summary, either the contribution of the off-shell gluon at the hard 
vertex is maintained, or it is compensated by the non perturbative input. 
In the first case (model 1), the hard photon is assumed to probe the gluon 
content of the target, and consequently, the 1/7^ singularity of the coeffi- 
cient functions is preserved. In the latter case (model 2), the photon probes 
the quark distribution and thus the coefficient function double pole can be 
compensated by the proton scaling function ^7(7) oc 7 when 7 goes to 0, see 
formula (12). 

As we shall see, a prototype of model 2 is provided by the wave function 



formulation of the photon-proton interaction [0, |I2[. In this framework, 
deep-inelastic scattering processes are formulated in terms of the probability 
distribution of a QQ pair (considered as a dipole configuration) in the vir- 
tual photon, convoluted by the dipole-proton cross section. In our case, the 
dipole-proton cross-section is described by the convolution of the probability 
distribution of primordial dipoles in the proton times the dipole-dipole BFKL 
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cross-section 0. This will allow a direct comparison between model 1 and 
2, which have a similar parametrization differing only by the pole structure 
at 7 = 0. 

In references [|ll], |T^, one finds the expressions for the wave function and 
probability distribution of the photon QQ states. The virtual photon can be 
described in terms of probability distributions (when the interference terms 
[ |T9[] are not relevant) depending on the quark mass M and charge e 

$;(., r; Q^ M') = {{z' + (1 - .)^)6^ir?(er) + M^Ko^(er)) , 

<!>l{z,r;Q^M') = ^^Ae'Q\\l-zfKl{er) , (22) 

where = z{l — z)Q'^ + M^, and the i^o,i are the Bessel functions of second 
kind |15|. $j> j;^(z, r; Q^, M^) are the probability distributions of finding a 
dipole configuration of transverse size r at a given z, the variable z (resp. (1 — 
z)) being the photon light-cone momentum fraction carried by the antiquark 
(resp. quark). 

The transverse and longitudinal total cross sections ct^.l read 
(Tt,l = I d^rdz ^(r, z; M^) f d\dzp $P(rp, Zp) ad{r, r^; Y) , 

(23) 

where we have introduced the probability distributions $^(rp, Zp) of dipoles 
inside the proton ||^, ^ The dipole-dipole cross section ad{r,rp]Y) is 
assumed not to depend on z. This hypothesis means that we neglect sub- 
asymptotic effects related to the momentum carried by the quarks (while we 
do not neglect the quark masses). In practice, it corresponds to consider 
iV ~ in the Mellin-transform with respect to rapidity as in section 2. 
adir,rp;Y) reads 0,^0: 



a,{r,rp;Y) = Anrrp J ^ (^) e"^>^W^^,(7) , (24) 

where x{l) is the BFKL kernel (p!OD and the elementary two-gluon exchange 
amplitude is given by 
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We define the non-perturbative scale Qo characterizing the average dipole 
size by ^0 



/dvr/^^p*^(-p,-p)-^^, (26) 

where Ueff can be interpreted as the 7-dependent average number of primary 
dipoles in the proton. Finally, the Mellin-transform of the photon wave- 
function is defined by: 

J^{r'y-'Jdz<l>l,{r,z-M') = <PrAr,M'){Qy-' , (27) 

where we have exphcitely factorized the photon scale dependence. After 
plugging these formulae into equation (|23|) , and performing the integrations 
with respect to r, Vp and z, Zp one finds: 



crT,L = / 77- I 7^ I e ^ ^(^)^0T,L(7;M2)A,K7)^e//(7) • 



Q2 J 2m \Q'o 



(28) 



Inserting the virtual photon probability distribution j;^(r, z; Q^, M^) into 
equation (|27D and after a straightforward but tedious calculation, one finds 
expressions for $t,l(7; M'^) which can be cast into the following form: 



0,,(7;M^) = ^^^^^^i^^^ /2-7+3(,_^).r(l-7) 
as 4:71 7 [ r(7) 



(29) 



where hT,L{'l'i M"^) H is related to the Mellin transform of a accounted for 
in the preceding section (see equation (3)). It is clear from formula (^) 
that the double pole of /iT,L(7; M"^)/"! is turned into a single pole due to the 
factor {...}. This shows that the factorization of the probability distributions 
4>t,l{1', M^) at the hard vertex leads to the singularity structure (^) of model 
2. 

Some comments are in order: 

• The 7-dependent but M^-independent factor {...} is nothing but the 
coupling of the virtual gluon to the QQ pair configuration of the vir- 
tual photon wave function. Indeed, the two-gluon exchange elementary 
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amplitude p5| ) can be rewritten as 



»2 

Ml) = = ^M.Hi-,) , (30) 

with 

2-2^-1 r(i-7) , , 

which is, up to a factor the eikonal couphng of a gluon to a dipole 



IP, ^ at the lower vertex. Thus, the factorized factor {...} is nothing 
else than f (1— 7), i.e. the coupling of the gluon to the photon's dipole 
configuration, as shown in figure 1. 



Using (|30| , |3TD , equation ( poD can be rewritten 



Aeiil) ■ (PtAt, M ) = aem— — f (t) • (32) 

In the hard perturbative domain where we consider a photon-dipole 
interaction with a dipole of small size, formula (p2D means that both 
high energy factorization and the wave-function formalism are iden- 
tical. The cross section can be equivalently factorized in two ways: 
either by the convolution of the photon gluon cross section times the 
gluon coupling to the dipole (right hand side), or by the probability 
distribution of a pair of quarks in the photon times the dipole-dipole 
elementary interaction (left hand side). 

However, as previously discussed, the non perturbative input may lead 
to distinguishable models. 

At the level of the non perturbative input, we note a relation between 
the two formulations (12) and (^), namely 

^(7) = ^e//(7) , (33) 

TT 7 

which generalizes the result for M = 0, obtained in ref. M. 
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5 Phenomenology 



Following our theoretical discussions of the previous sections, we will consider 
two definite models relying on different formulations of the residue function 
c<j(7). On the one hand, the model 1, with 

uj{j) = Ci (constant) (34) 

corresponds to the factorization at the gluon level (see figure 1). On the 
other hand, in the wave function formulation of the dipole model of section 
4, we can reformulate the integrand of the structure function (see equation 
(12)) by using relation (|32[), namely: 

^-(7)= (;^^) ^(7) -^(7) -.(1-7). (35) 

In the model 2, the hard vertex is thus described by factorizing ^(7) which 
means that we consider 

c<j(7) ■ t'(l— 7) = C2 (constant) . (36) 

As expected, oo{'-)) behaves like 7 when 7 goes to zero in this framework. 
In the large I7I region (7 = 1/2 + iz/, v goes to infinity), one has |li^(7)| — 
|C2/f(l— 7)1 oc |7p. The exponential cutoff of 0r,L is the same as the one 
(19) of the coefficient function and thus the large I7I constraint is satisfied. 
Note that, in both cases, the proton structure functions depend on three free 
parameters only: the global normalization Ci or C2, the effective constant 
strong coupling a^, and the non perturbative scale Qo- 

We determine these parameters for the two models by a fit of F2 = Ft+Fl 
in their kinematical region of validity {x < 10~^). In this region, is 
automatically limited by the HERA kinematics. Using the corresponding 
103 experimental points given by the HI collaboration ||2^, we fit our results 
with the contribution of the three light quarks u, d, s (assumed massless) 
and of the charm quark (mass Mc). As usual, we will vary in the range 
1.35- 1.7 GeV §. 

The F2-fit for the medium mass Mc = 1.5 GeV is displayed in figure 2, 
together with the predictions for its charm component F2. In table I, we give 
the cind the value of the fit parameters for the fits for model 1 and 2 and 
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for Mc = 1.35, 1.5, and 1.7. For model 1, the per point is always less 
than 0.9, while for model 2 it is even lower. 

Some comments on the parameter values are in order. For fit 1, the value 
of Qq is around 330 MeV which is a typical non perturbative scale for the 
proton. The value of the effective coupling constant in the BFKL mechanism 
as (0.07) is rather low. This value would amount to an effective pomeron 
intercept ap = 1 + asA^c41og2/7r ~ 1.18. (It is known that ap is influenced 
by sub-leading corrections to the BFKL kernel, see e.g. references [^). Note 
that the fit of model 1 for F2 in the framework of the QCD dipole model and 
high energy factorization is compatible with the previous ones in the same 
framework @. 

We also in parallel performed the phenomenological analysis using the 
model 2-ansatz (|36D. The obtained fit reproduces fairly well the data for F2 
(see figure 2). As indicated by the lower value of x^, it seems that the data, 
especially in the small-x small-Q^ region are even better reproduced than for 
model 1. Note that the value of Qo — 1.2 is substantially higher and the 
effective coupling constant as is larger (~ as{Mz))- 

Following the factorization scheme and using formulae (^ (model 1) or 
( p3D (model 2), we obtain now parameter- free predictions for the longitudi- 
nal structure functions Fl and F£ (see figure 3). The predictions for Fl and 
F£ are shown together with the indirect HI determination of Fl [^. As 
expected, F£ becomes a significant part of Fl at small x and high Q^. As 
already noticed P], the predictions are low but compatible with the present 
large error bars. Note also that the determination depends on the theo- 
retical scheme one considers p6|, so it is difficult to draw any conclusion on 
Fl at this stage. 

The main outcome of our analysis is a parameter free prediction for F|, 
the charm component of the structure function (see fig.4). When compared 
to ZEUS and HI data for the charm component in the same range in x and 

0, we find a good agreement within the present experimental error bars. 
The extrapolation of the prediction to the kinematical range of EMC data is 
correct while the dipole model is not expected to be valid in this region. 

Looking in more detail to the predictions of model 1 (figure 4-a) and 
model 2 (figure 4-b), we observe the following features. The dispersion of 
the results with respect to is rather small for fit 1 with a maximum of 
10% when x reaches 10~^. Moreover, the prediction is comparable to the 
next-leading order GRV prediction which proves that F2 cannot allow one to 
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distinguish between the two approaches. 

For model 2, the predictions obtained for F| are displayed in figure 4-b. 
The prediction for the HERA region is quite satisfactory and, interestingly 
enough, somewhat higher than both the NLO GRV predictions and fit 1. 
It would certainly be useful to deserve some experimental and theoretical 
attention to this difference. Indeed, the scheme of model 2 which is based on 
a different factorization than both GRV and model 1 could lead to a better 
understanding of the data. A sensible decrease of the experimental errors on 
F2 in the region where is moderate (of the order 10 GeV^) and x small 
(of the order 10"'^) could allow a refinement of this analysis. 

All in all, the comparison of fits 1 and 2 shows that the high energy 
factorization prediction for F2 depends on the non-perturbative input in 
the HERA range. Indeed, a different factorization scheme such as model 2 
corresponds to a modification of the non perturbative input of high energy 
factorization. However, the largest uncertainty due to this effect is of the 
order 20% and less than the experimental uncertainty. Thus, there is no 
present evidence of a distinction between these schemes and the DGLAP 
scheme. 

6 Conclusions and outlook. 

Let us summarize the results of our analysis. 

(i) High energy factorization [|1| gives predictions for the x and behaviour 
of F2 at low X in the framework of models based on the BFKL dynamics. The 
charm prediction is fixed by high energy factorization once the total proton 
structure function F2 is fit. More generally, this prediction is well defined for 
any quark mass. 

(ii) The factorization of the hard virtual photon vertex at the level of the 
exchanged gluon (model 1) or the exchanged quark (model 2) leads to two 
specific classes of models which both give a satisfactory fit of F2. The for- 
mulation |TI|, based on the virtual photon wave function is shown to lead 
to a physical realization of model 2. 

(iii) Both models lead to predictions for F2 in agreement with the present 
data. However, the second scheme leads to a higher F2 than both model 1 
and NLO DGLAP predictions P]. This justifies further experimental and 
theoretical studies. 
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(iv) The derivation leads to a new constraint on the unintegrated gluon dis- 
tribution in the proton due to the renormahzation group evolution at high 
transverse momentum of the off-shell gluon. 



It is interesting to discuss the gluon distribution functions in this frame- 
work. In the high energy factorization scheme, it would consist in an ex- 
pression like (12), replacing the coefficient functions by unity. However, the 
convergence of the integral (12) is then not a priori preserved and even if 
convergent, the dependence of the result on the large I7I region would be 
larger, casting a doubt on the relevance of the BFKL dynamics. In the wave- 
function framework, the extraction of the gluon structure function is even 
more problematic since the factorization is at the quark level (see fig.l). On 
a more general ground, this confirms the statement of caution [jl^ about 
extracting the gluon from this kind of analysis, especially at low x. We think 
that this point deserves more studies. 

Finally, the running of the coupling constant and other aspects of next 
leading order BFKL resummation features have been neglected in the present 
analysis. The small value obtained for the effective coupling constant in the 
fits clearly indicates that such effects are important and should be included 
in a more detailed analysis. Indeed, a preliminary theoretical hint on the 
effective behaviour of the pomeron singularity at NLO accuracy leads to a 
small and constant intercept of the order 0.2 A NLO BFKL analysis 



of F2 would thus be required in view of a future improvement of the experi- 
mental analysis. 
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TABLE CAPTION 



Table I 

Total proton structure function: ^"'^ parameters for the fits 

For each model, we give the value for 103 points and the parameters 
C, Qo and as- Three different charm masses, Mc — 1.35, 1.5, 1.7 GeV are 
considered. I-a. Model 1. I-b. Model 2. 
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Ci 


Qo 


as 


1.35 


77 


19.88 


0.334 


0.070 


1.50 


82 


20.27 


0.326 


0.070 


1.70 


88 


20.70 


0.321 


0.070 



a. Model 1 







C2 


Qo 


as 


1.35 


51 


8.68 


1.208 


0.105 


1.50 


51 


9.22 


1.215 


0.104 


1.70 


52 


9.86 


1.227 


0.103 



b. Model 2 



Table I 
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FIGURE CAPTIONS 



Figure 1 Quark- antiquark pair leptoproduction at high energy 

The upper vertex represents the virtual photon (with momentum q)- 
virtual gluon (with momentum k) fusion diagram through the production 
of a pair of quarks of mass M. The proton (with momentum p) interacts 
with the gluon through the BFKL kernel. The non perturbative proton ver- 
tex is schematized by the shaded area. Model 1 corresponds to factorization 
at the gluon level, model 2 at the quark level, see text. 

Figure 2 The fits (model 1 and model 2) for the structure function F2 

The structure function F2 and its parametrization are displayed for model 
1 and model 2. The fits have been performed with the 1994 HI data (tri- 
angles) but only using the 103 experimental points for which x < 10~^. 
Continuous line: fit for model 1; dashed line: fit for model 2; dotted line: 
prediction for the charm component F2, for Mc = 1.5 GeV and model 1. The 
available HI data are marked by stars. 

Figure 3 Predictions for the longitudinal structure functions ( continuous 
line) and F£ (dashed line) 

Continuous line: model 1; dashed line: model 2; dotted line: F£ (model 
1). The experimental points available from HI have been reported on 
the graph. 

Figure 4 Predictions for F2 

The reported data on the plots are from HI (squares), HI D* (circles), 
ZEUS D* (full circles), and at lower energy, EMC data (crosses). The predic- 
tions of our models are displayed by a band dehmited by the two continuous 
lines (Mc = 1.35 for the higher curves and Mc = 1.7 for the lower curves). 
Extrapolations of our predictions down to the EMC range are plotted as 
dot-dashed lines. 

Figure 4-a: model 1 predictions. The GRV prediction |^ based on NLO 
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analytic calculation is indicated by the shaded band (borrowed from ref. 
second paper). 

Figure 4-b: model 2 predictions: the band of solutions is delimited by thick 
lines. The model 1 predictions are recalled by the band delimited by thin 
lines, for comparison. 
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Figure 2 
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Figure 3 
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